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ABSTRACT 

The purpose of this thesis is to provide the student 
of queuing systems with a vehicle through which a better 
understanding of the interrelationships between stochastic 
processes and queuing systems can be achieved. 
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I. INTRODUCTION 



A queuing system is comprised of action facilities 
(servers) and waiting facilities (queues or lines) . The 
relationship between the components can be of many forms, 
each with their own characteristics. Consider a tool 
crib in a machine shop where the machinists are obliged 
to check out their tools; or consider the telephone system, 
in which callers request communication capacity. Both 
situations involve requests for action from limited resources, 
in these cases, servers or telephone circuits. In order to 
allocate resources properly, it is necessary to understand 
the interaction between requests for action and the actions 
themselves. Individually these can be modeled as stochastic 
processes, but in order to predict flow through the queuing 
system, it is necessary to understand the interrelation 
between the patterns of event occurrences. When demands 
occur at regular intervals and actions are of fixed dura- 
tion, it is relatively easy to predict flow; if there are 
more requests per unit of time than actions over the same 
time span, the queue will grow without bound unless the 
event relationship is changed. When variation is added to 
the timing of events, the flow is not as obvious. For 
example, if the average number of requests for action 
equals the average number of actions, per unit of time, but 
both occur with variability, then it would be useful to be 
able to predict such quantities as expected delay and action 
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facility idle time. Queuing models are used for this 
purpose . 

The purpose of this thesis is to provide the student 
of queuing systems with a vehicle through which a better 
understanding of the interrelationships between stochastic 
processes and queuing systems can be achieved. 

This is done through a combination of a computer simula- 
tion and a computer-driven graphical display. 
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II. QUEUING MODEL SIMULATION 



This computer program simulates queuing systems in two 
stages. First, a next-event type computer simulation creates 
a time history of the queuing system. Second, a dynamic 
graphical display depicts the occurrence of the queuing 
process in real time. 

Five basic, academically interesting, queuing models, 
each with parameter options, can be simulated. Queuing 
models are generally described in terms of the distribution 
of request occurrences (arrival distribution) , the distribu- 
tion of action occurrences (service distribution) , and the 
number of action facilities (servers). Specifically, the 
times between request occurrences are taken to be indepen- 
dent and identically distributed realizations from the 
arrival distribution. The duration of each action occur- 
rence is also taken to be an independent and identically 
distributed realization from the service distribution. 
Variations on these assumptions occur for some models and 
are described in Section A. 

The abbreviation G/G/C describes a queuing system in 
which the arrival and service distributions are general 
(unspecified) and there are C servers. 

The five options in this simulation are; 

(1) G/G/C 

(2) G/G/C with losses 

(3) G/G/C with feedback 
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(4) G/G/C with finite source 

(5) Two simultaneous G/G/C queues for comparison. 

The specific choices of arrival distribution and service 
distribution which can be simulated are fully described 
in Section B. 



A. THE QUEUING MODELS 
1. G/G/C 




Picture 1 



This is the simplest of the queuing models considered 
here. Each arrival waits in the queue until a server is 
available, enters service for a calculated period of time, 
then departs the system. 

This model may be simulated with from one to six 

servers . 
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2. G/G/C with Losses 




Picture 2 



Consider making a long distance telephone call on 
Mother's Day. I£ all the trunk lines are busy, your request 
for service will be rejected. The G/G/C model with losses 
reflects this by rejecting arrivals who find the queue 
full. The queue limit used is specified as an input 
parameter. Note that rejected arrivals are lost and are 
not assumed to return later. The fraction of all arrivals 
who are served (not rejected) is recorded for subsequent 
display . 
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3. G/G/C with Finite Source 




Picture 3 



Consider a closed-circuit educational television 
system at a school. The finite population consists of all 
the televisions in the system. When one fails, it enters 
the queue to be serviced. After repair, it re-enters the 
population pending the next failure. 

This situation is modeled by assuming that each 
action source, upon entering the idle population, is 
assigned a waiting time from the arrival distribution. 

When this time expires, the source requests service, 
joining the queue if necessary. As each action source 
terminates its service, it returns to the idle population. 
The user specifies the number of servers (one to six) and 
the size of the population. 
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4 . G/G/C with Feedback to the Queue 




Picture 4 



This is an extension o£ the G/G/C model, with from 
one to six servers, allowing customers leaving service to 
depart the system or to return to the queue. At each 
service completion, the customer leaves the system with 
probability P or returns to the queue with probability 1-P. 
The user specifies P as an input option. 
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5 . Two Simultaneous G/G/C Queues 




Picture 5 



This option allows comparison o£ two standard 
G/G/C queues, each with up to three servers. All parameters 
o£ each system are speci£ied separately, and, i£ desired, 
identical arrival streams may be utilized. This allows 
comparisons between di££erent arrival and/or service distri- 
butions and associated parameters, and between a di££erent 
number o£ servers. Also, the e££ect o£ the various queuing 
disciplines can be compared. 

B. ARRIVAL/SERVICE DISTRIBUTIONS AND QUEUE DISCIPLINES 
As well as options already described, the £ollowing 
arrival and service distributions are available to the 
user o£ the program. Any o£ the £ollowing probability 
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distributions can be used as an arrival or service 
distribution. 

1 . Distributions 

a. K - ERLANG (exponential when K=l) 



F(t) 



C3t)^ g-3t 



E 

J=0 



J! 



RATE X = I 

CV = — 

/K 

Input parameters are: 

(1) K 

i 

(2) Rate - in terms of mean number of events 

per minute. 

b. Hyperexponential 

F(t) = p(l-e'^l^) + q(l-e’^2^) 
where : 

(1) p + q = 1 

(2) CV = — 

/p 

(3) X = Rate is specified and then 
X^, chosen to satisfy-- 

X 2 = (2 + /I)X 
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I 




X 



1 



(PC2_L_r^KL) 

2 + /2 - q ^ 



Input parameters are: 

(1) Rate, specified in terms of mean number 
of events per minute. 

(2) Coefficient of variation (must be greater 
than 1.0) 

c. Degenerate (Fixed Rate) 



0 t < 1/X 

F(t) = 

1 t > 1/X 

Input parameter: 

(1) Rate X specified in terms of mean number 
of events per minute. 

2 . Queuing Disciplines 

For comparison purposes, there are four queuing 
disciplines to choose from. These determine the priority 
by which members of the queue will enter service. 

a. FIFO - the first into queue is the first out. 

Over all disciplines which are independent of service times, 
this one minimizes the variance of delay. 

b. LIFO - the last into the queue is the first 

out. Some systems, by their nature, function in this manner. 

c. SSTF - the queue member with shortest service 
time is served first; this minimizes the mean delay over 
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all other disciplines, but the variance of delay is greater 
than that of FIFO, 

d. RANDOM DRAW - applies to a system in which 
insufficient information is known to implement any other 
discipline . 

3 . Length of the Simulation 

The period of simulation is specified in minutes 
by the user. All queues are empty and servers idle at the 
start of the simulation. 
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III. THE METHOD 



A. THE COMPUTER SIMULATION 

The computer simulation is a modified next event type. 
There are two types of events, arrivals to the queue and a 
completion of service. With the first there will be an 
associated entry to service if there is an idle server, 
and with the last an associated entry to service if the 
queue is not empty. 

The simulation is begun with the generation of the first 
arrival time and its associated service time. The times are 
generated via the transformation of random numbers to stochas- 
tic variates of the type specified. At the occurrence of 
an arrival, the next arrival time and associated service 
time are generated. When a server is available, a waiting 
member of the queue enters service with completion time 
equal to the sum of his assigned service time and the time 
he departs the queue. Thus with the clock starting at zero, 
a time history of events is created which tells what happened 
and when. 

At the same time, the appropriate values are retained 
from which statistics of the simulation can later be 
calculated . 

B. THE GRAPHICAL DISPLAY 

At the completion of the simulation phase, the time 
history is depicted dynamically on the cathode ray tube 
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in real time. Pictures 1 to 5 show the form of the display, 
with the man-like images in motion along the appropriate 
path, and the number in the queue box showing the queue’s 
current length. 

The timing and motion of the display are controlled 
by the analog computer via the hybrid connections to the 
digital computer. The first connection is for the purpose 
of timing until the next event is to be initiated. In this 
case, the analog computer integrates from zero to a preset 
time value in volts, at the rate of one volt per second. 

The interrupt to the appropriate subroutine then occurs, 
and the analog computer is reset with the next time value, 
after which the integration is reinitiated. 

The second connection causes the dynamic motion to 
occur by continually checking events in progress and causing 
the position of the corresponding images to be incremented 
appropriately; this interrogation and incrementation occurs 
approximately 10 times per second depending on the settings 
of the analog computer. Movement increments are between 
one-half inch and two inches (depending on the path) , 
causing a motion of between five inches and 20 inches per 
second. 
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IV. STATISTICS FROM THE SIMULATION 
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2. 


PARAMETER LISTING 



At the end of the simulation period, the following 
statistics are calculated and displayed, as shown above 
AT = time - average length of the queue 

BT = time - average number of busy servers 

LT = time - average number in the system = QT + BT 



QC = customer - average length of the queue 

BC = customer - average number of busy servers 

LC = customer - average number in the system = QC + BC 

DC = customer - average delay in the queue 

SC = customer - average service time 

WC = customer - average wait in the system = DC + SC 
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DM = maximum delay experienced by any customer 
DO = probability that a customer's delay is greater than 
zero , 

DX = probability that a customer's delay is greater than 
XD. (The value of XD is an input option.) 

QM = maximum queue length reached during the simulation 
QO = probability that the queue length will be greater 
than zero. 

QL = probability that queue length will be greater than LQ 
(the value of LQ is an input option) . 

DV = standard deviation of the delay in queue 

SV = standard deviation of the service times 

QV = standard deviation of the queue length 

NC = total number of arrivals 

PS = percentage not lost (losses model only) 
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V. PARAMETER INPUT 





I TYPE 
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RUNTIME 


1 .tt 
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l.tt 


1 .ff 
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l.tt 


lOSPuN 
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1 


HAXO 
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XO 


1 .ff 


NPOP 


1 




I DUAL 


1 


P 


l it 



ENTER 1 . CONTINUE 

2. STATISTICS 



Picture 7 

Parameters can be input by the card reader or by the 
teletype terminal. In either case, the namelist form is 
used, where the entry form is NAME = X, NAME being the 
variable title and X being the value being assigned. For 
real variables, the decimal point must be entered. Name- 
list input is terminated with a *. 



A. SAMPLE INPUT CARD 
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B. NAMELIST INPUT VIA THE GRAPHICS TERMINAL 

If this mode of entry is specified, one variable at 
a time is entered, followed by a carriage return. Entries 
are displayed, as entered, on the graphics screen. Termina- 
tion occurs with the entry of a * followed by a carriage 
return . 

C. INPUT VARIABLES 

The following conventions are employed for specifying 
the options and parameters to be used. 

ITYPE = 1 gives G/G/C 

2 gives G/G/C with losses 

3 gives G/G/C with finite source 

4 gives G/G/C with feedback 

5 gives two simultaneous G/G/C queues 

RUNTIME = desired length of simulation (minutes) 

IC = number of servers for queue 1 

ICl = number of servers for queue 2 (only applicable if 



ITYPE 


= 5 is chosen) 


IDIST(I) = 1. 


K - ERLANG 


2. 


HYPEREXPONENTIAL 


3. 


DEGENERATE 



KK(I) = K in K-ERLANG (use K=1 for exponential) 

RATE (I) = Rate in number per minute 

CV(I) = Coefficient of variation (specify for IDIST(I) = 2 
only. This variable is computed otherwise.) 
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where I 



= 1 for arrivals to queue 1 (except finite model) 

2 for arrivals to queue 2 or finite population 

idle parameters 

3 for service for queue 1 (except finite model) 

4 for service for queue 2 or finite population 

service parameters 



IDSPLN(I) 



where I = 



= 1. FIFO 

2. LIFO 

3. Shortest service time first 

4. Random draw 

1 for queue 1 

2 for queue 2 



LQ = parameter in the statistic P (Q > LQ) 

XD = parameter in the statistic P(wait > XD) 

IDUAL = 0. Independent arrival streams 
1. Identical arrival streams 
(applicable to ITYPE = 5 only) 

MAXQ = maximum queue size (for G/G/C with losses only) 
NPOP = population size (for finite source option only) 
P = completion probability (for feedback option only) 
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VI . THE PROGRAM 



The program requires the use o£ the XDS 9300 digital 
computer, the Cl 5000 analog computer, and a graphics 
terminal using the GATED monitor. 

The analog computer wiring diagram necessary to drive 
the display sequence appears in APPENDIX A. Before 
beginning, set Delay Flop Switches DFOO, DFOl , and DF02 
to . 1 MS and set counter FOO to 1 and F04 to 0. Then select 
RUN and DIGITAL COMPUTER on the analog control panel. This 
provides for movement interrupts to occur at a rate of 
100/11 hertz. If desired, the timing values will be dis- 
played on the RATIOMETER by selecting TRUNK 420 prior to 
DIGITAL COMPUTER. 

The program consists of a FORTRAN program and a META- 
SYMBOL 9300 program. Because the history is buffered out 
to the device assigned the logical unit value 7, an addi- 
tional control card is needed. 

Job cards for deck run should be: 

"I JOB 
-lAGT 

tASSIGN 7=DF1A 
-iFORTRAN LS,GO 

FORTRAN' DECK _ 

-iMETA 9300 SI,LO,GO 
META- SYMBOL DECK 
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-lEOF 



-iLOAD XR,MAP 
-I DATA 
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ffpro 










PROGRAM STEPS 
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APPENDIX A 



ANALOG SCHEMATIC 



SCOPE Y2 





T 420 
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CLOCK( I3)=TIME-PSTI ME 
KIM( I3)=N0(2) 
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45 CALL TIMES 

41 ET IME{ I+6)=TIM(KP( I ) + l, I ) 

42 IF(LINE( I) .GT.MAXLNd ) ) MAX LN( I ) =L IN E ( I ) 
RETURN 
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